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I. INTRODUCTION 



Scintillation is the term used to describe the 
atmospherically induced intensity fluctuations of a constant 
power output laser source as seen from a distant detection 
device. The precise manner in which a detector’s aperture 
size affects the measured value of the variance of the 
intensity fluctuations as a function of distance from the 
source and lasing wavelength has not been satisfactorily 
established. Several aperture averaging experiments have 
been conducted under various turbulence conditions to examine 
this effect along with an analysis of the frequency power 
spectrum effects of aperture size. 

A. PROJECT HISTORY 

Since 1971 an interdisciplinary group of professors 
principally from the Departments of Physics and Chemistry, 
Meteorology, and Mechanical Engineering have been investi- 
gating under a Navy research contract various aspects of 
laser beam propagation over the marine boundary layer. Data 
has been gathered thus far on (1) Modulation Transfer Function 
(MTF) at wavelengths of .4880 ym, .6328 pm, 1.06 ym and 10.6 ym 
(2) scintillation at .4880 ym, .6328 ym and 10.6 ym, (3) extinc 
tion at .6328 ym, (4) beam wander and (5) beam spread. 

The stated general objectives of the project are: (1) to 

develop an understanding of the factors affecting the propaga- 
tion of optical radiation through the atmosphere, particularly 
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the marine boundary surface layer, the lowest 30 meters 
above the ocean surface; (2) to develop the capability of 
predicting optical system performance on the basis of point 
meteorological measurements. 

Previous thought had been that the sea surface might 
present a smooth, isothermal layer, giving rise to little 
turbulence effect on a propagating light beam. However, 
early experiments in the project proved that this was not 
the case. 

B. SCINTILLATION 

1 . Description 

The twinkling one sees when observing distant light 
sources is a physical realization of scintillation. This 
appears as fluctuation in the output from a small photo- 
detector illuminated by a distant source. 

In the general case of finite beam propagation, 
measured intensity fluctuation can be caused by any or all 
of the following: (1) source output power fluctuations, 

(2) beam wander, (3) detector motion through a plane 
perpendicular to the optical axis, (4) a temporal wavefront 
breakup within the beam caused by small-scale spatial and 
temporal variations of the dielectric constant (or index of 
refraction) of the lower atmosphere. Of these, the most 
critical for experimentation or information transmission are 
the turbulence effects of scintillation, beam wander and 
beam spread, since they cannot be overcome by system 
refinement . 
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2. The Effects of Turbulence on a Propagating Wavefront: 
Physical and Theoretical Models 

a. Kolmogorov Model of Turbulence 

Currently the only well-developed model of 
atmospheric turbulence is that of Kolmogorov. It is basically 
a velocity turbulence model that assumes homogeneous, iso- 
tropic turbulence. The turbulence is assumed to be composed 
of small turbulent eddies, or turbules, each of characteris- 
tic temperature and density (hence index of refraction) , as 
illustrated in Figure 1. In this homogeneous turbulence 
model, often referred to as "frozen-in turbulence," fluctua- 
tions in the dielectric constant over a vanishingly small 
volume element are due to atmospheric temperature variations 
produced by the turbules being swept past the reference 
point by the mean wind. 

Let n(r, t) be the refractive index at point r 
in space at time t. Then 



where n^Cr^t',!) is the deviation from the average value at 
point r and time t. 

The refractive index structure function gives 
the mean square variation in the refractive index between 



x 



CD 



n(r,t) 




o 



= <n(r,t)> T + n 1 (r,t ’ ,x) 
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two points in space: 



(2) 


D n (p ) ■ < [ n ldl) • n id 2 )] 2> . 


where 


(XJ 

f-*| 

1 

mT 

II 

Ql 




The cross -correlation function of refractive 



index fluctuations is defined as 



(3) 


T 

B n (p) = \ J n 1 (r 1 ,t')n 1 (r 2 ,t')dt' 
0 




= <n 1 Cl 1 ) n iCr 2 )> 



Under the assumption that homogeneous turbulence is present 
and the appropriate statistical averages for B n (p) and 



D n (p) are 


stationary in time and locally constant in space, 


one finds 


that 


(4a) 


D n (p) = 2[B n (0) - B n (p) ] 


(4b) 


B n (P) - £[D n C») - D n (p)] . 



As p becomes smaller, correlation increases while the 



structure 


function decreases. 



According to the Kolmogorov theory of turbulence, 



for £ << 

0 

scales , 


p << L , that is, within the "Inertial Subrange" of 
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( 5 ) 



) 




and the corresponding temperature structure function for 
meteorological measurements is 



& Q and L q are the inner and outer scales of turbulence 
in the Inertial Subrange (Figure 2) , and and C,p are the 
refractive index structure constant and temperature struc- 
ture constant, respectively. If l (or L) represents the 
diameter of a turbulent eddy, i is the scale size below 
which viscous effects cause the conversion of turbulent 
energy into heat, and subsequent destruction of the eddy 
ensues. L , the outer scale size, is the minimum scale 
at which energy enters the turbulent system from wind shear 
and heat flux. 



(the von Karman form) of the power spectrum (the three- 

dimensional Fourier transform of the cross -correlation 

function) , which illustrates its approximate dependence on 

l and L , is 
0 o 



( 6 ) 



d t (p) = C T 2 p 2/3 



For Kolmogorov turbulence, an analytical form 



(7) $ n (K) = 0.033 C R (K + L q ") 



exp [ - (K& q ) 2 ] . 
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For K >> L q , this reduces to the more common form used by 
Tatarski , 



( 8 ) 



MK) = 0-033 C„V 11/3 exp[-K 2 /K m 2 ] , 



n 



n 



where 



v _ 2tt 

K m z 

o 



b. Physical Models 

(1) Tatarski 's Geometrical Optics Approach . The 
physical mechanism through which turbulence induces beam 
breakup is by the focusing of sections of a wavefront 
produced by the inhomogeneities in the index of refraction 
over the propagation path, Tatarski uses a plane wave 
incident upon a series of closely packed refractive spheres 
of various sizes. Using the paraxial approximation, a 
sphere acts as a lens with focal length f = , where n^ 

is the index of refraction perturbation within the sphere. 



For a = 1mm, Z Q of the Inertial Subrange, and n^ = 10 



-6 



f = 500m. For a = 1m, the lowest L limit, f = 500 Km. 

Thus, for most applications in a turbulent medium, a detector 
will appear to be closer to the "lens" than to the focal 
plane. The fractional intensity change induced per sphere 
will be small. Over a long path, perturbation effects will 
be multiplicative and the resulting wavefront at the 
receiver will be neither coherent nor of constant intensity. 
In addition, diffraction effects, important at distances 
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2 

greater than the Rayleigh range, Z R = ka , limit the ability 
of the sphere to focus the incident wave. The amplitude 
fluctuation induced by a sphere at a distance Z from a 
detector is estimated [Ref. 1] as 



(9) 



6A ~ Z _ 



2 I n jJ Z 
a 



, Z << f . 



From equation [9) it is seen that the largest fluctuations 
are produced by the smallest turbules. By the geometrical 
optics approach [Ref. 1], within the Rayleigh range, 
amplitude fluctuations increase as the cube of the distance, 
are dominated by the smallest turbules, and are independent 
of wavelength: 

CIO) <(|V> ' C N 2 * " 7/3 z 3 . 

O 

However, this model does not take into account diffraction 

effects, important for smaller turbules that contribute most 

toward amplitude fluctuations. 

(2) Diffraction Effects . Diffraction effects 

are taken into account by assuming that spheres of radius 

~ (XZ») ' , where Z' is the distance from the center of 

the turbule to the receiver, contribute very weak, diffracted 

waves with no focusing effects at the receiver. Substituting 
1/2 

( A Z ) for in equation (10) , the mean square amplitude 

fluctuations become 
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2 k 7/6 Z 11/6 , kl 2 
* o 



< Z < kL 



o 



2 



6A 

Now defining where X-^ is the first order log- 



amplitude for | Xj | << 1, Tatarski has calculated the first 
order log-amplitude variance with diffraction effects con- 
sidered for plane waves as 



state that one would expect the behavior of a beam wave to 
lie between that of the plane and spherical waves. There- 
fore, one may use equations (12) and (13) as bounds for 
weak scintillation. For inhomogeneous turbulence, equations 
(12) and (13) can be expressed in the form 



( 12 ) 




and for spherical waves as 



(13) 




2 

both restricted to Z >> k& Q . Lutomirski and Yura [Ref. 1] 



Z 



(14) < Xl 2 > (P) = 0.56 k 7/6 




o 



Z 




C N 2 (Z')dZ* , 



o 
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where Z is the distance from source to receiver and Z' 
from turbule to receiver. For plane wave propagation, the 
fluctuation effects are weighted by the factor (Z-Z*)^^, 
which has its maximum value near the source and approaches 
zero at the receiver. The spherical wave weighting factor, 

[ (Z-Z * ) (y—) ] , has its maximum value at the midpoint of 
the path and approaches zero at both ends . 

c. Theoretical Models, The Scalar Wave Equation 
Numerous attempts have been made to solve the 
scalar wave equation to relate the electric field E(r) 
of an optical wave passing through a turbulent medium to 
the refractive index fluctuations. Three approaches current- 
ly in use for solving this problem are: (1) Geometrical 

Optics, (2) The Born Approximation, and (3) The Rytov 
Approximation. Below is a brief resume of the approaches, 
their assumptions, highlights, methods of approach and 
restrictions. All three methods start with the scalar wave 
equation, 

(16) V 2 E + k 2 n 2 E = 0 . 

The derivation of this has assumed a medium with no free 
charges or current density, with the permeability of free 
space, that correlations between the components of the 
field polarized in different directions can be neglected, 
and that X << & . 
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(1) Geometrical Optics Approach 



(a) Amplitude fluctuations are attributed 
to the focusing and defocusing of the rays by the curvature 
of the turbulent eddies along the path. 

(b) Phase fluctuations are calculated from 
the changes in velocity as a ray passes through regions 

of different refractive index. 

(c) Diffraction effects are ignored. 

i/7 

(d) The major restriction is (XZ) << i . 

This reduces the applicability of this model to extremely 

1/2 

short ranges. Note: When (XZ) >_ Z , diffraction effects 

become important. 

(2) The Born Approximation 

(a) This is a perturbation technique to 
incorporate diffraction effects into the scalar wave equation, 
n = <n> + n-^ is substituted into the scalar wave equation, 
which is solved for the first order solution for E. The 
result is 

(17) E 1 = J 2k 2 ni (r')E o (r’)G(|r-r’ |)dv' , 

v 

where G ( | r - r 1 | ) is the free space Green’s function. 

(b) The perturbation technique used is valid 
as long as the amplitude fluctuations remain small. 

(c) It is assumed that the field at the 
receiver is the sum of the original incident field and the 
field scattered one time from a turbule. 
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1/7 

(d) It requires that (AZ) ' > . 

(e) Restrictions: In the optical wavelength 

range, amplitude fluctuations often become very large 

compared to the average value of the amplitude. Multiple 

scattering effects cannot be ignored (but are). References 

2, 3 and 4 contain details on multiple scattering effects. 

(3) The Rytov Approximation (Method of Smooth 
Perturbations J 

(a) This is a perturbation technique applied 
to a transformation of the scalar wave equation using the 
exponential notation: 

E (r) = e^ r ) , t|>(r) = X(r) + i<J>(r) . 

i(/(r), the "complex amplitude," contains both the amplitude 
and phase terms. Using this, the scalar wave equation takes 
the form of the non-linear Ricatti Equation, 

(18) V 2 \p + V\p-Vrp + k 2 (1 + 6n) = 0 . 

Now applying the perturbation technique 4* = l P Q + ty-p where 
is the perturbation to the average value of \p Q , one 
obtains the solution, 

(19) ^(r) = £-1^ J G(r-r') + k 2 6n]E Q (r ' ) dV 1 . 

0 v 
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The M Rytov Approximation" is = | Vip ^ | 2 = 0, and 

results in 

* k 2 ' r 

(20) * 10 (r) = J G(r-r«)6n(r»)E o (r’)dV» , 

0 v 

the "first order complex phase." A higher order solution 
for can be obtained by iteration. 

(b) This approach includes multiple scattering 
effects. Caveat- -these effects are included in an inflexible 
way which .does not depend on the turbulence or other corre- 
lated parameters [Refs. 5 and 6]. 

(c) The method requires that the perturbations 
of the log-amplitude, the real part of the complex phase 

A 

[X x = log e V 1, and of the phase, <J>^, be small, i.e., 

X|X 1 | << 1 and A|<J>^| << 1. 

2 

(d) Restrictions: A << , | Vi|> q | << |k Sn| , 

« !. 

I 

(e) This method includes both the Geometrical 
Optics approach and Born Approximation as special cases. 

(f) Theoretical evidence indicates the 

equation derived for amplitude fluctuations is valid as 
2 

long as <x > <_ 0.8 (Saturation Condition) [Refs. 7-10]. 

3. The Form of the Probability Distribution for 
Amplitude Fluctuations 

The application of the Central Limit Theorem in 
the derivation of the log-amplitude equation leads to a 
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prediction of a log-normal distribution of the amplitude 
fluctuations, or equivalently, a normal distribution for 
log-amplitude. Some authors [Refs. 11 and 12], taking a 
different approach, argue for a Rayleigh distribution. The 
Rayleigh distribution arises when the field at the receiver 
is the sum of the electric field components from randomly 
scattered wavefront elements. This occurs in line-of-sight 
propagation when the effects of the turbulent medium are 
far from the receiver. In line-of-sight propagation where 
the receiver is located within the region in which turbulent 
effects are being introduced, the field at the receiver is 
the result of multiplicative effects on a sin>gle wavefront. 
Strohbehn [Ref. 2] shows that i|> r (t), representing the ampli- 
tude and phase terms of a single wavefront, is composed of 
the sum of a large number of independent components, and 
from the Central Limit Theorem, has a normal or Gaussian 
distribution. 

For long enough paths, however, both the multiplica- 
tive effects on a single wavefront and additive effects 
from different wavefronts may be important, and Strohbehn 
indicates that it is difficult to tell which effect will be 
predominant. Empirical results thus far support a log- 
normal distribution [Refs. 14, 15 and 16 and Figure 9], 
However, extremely long-range experiments have not as yet 
been conducted. 



22 



C. THE SATURATION OF SCINTILLATION 
1 . Description 

For plane and spherical waves (only) in homogeneous 
turbulence and for weak fluctuations. 



( 21 ) 





* 4 <X 1 > 







2 . 2 
where is the normalized variance and a ^ the log- 

intensity variance. In terms of the log- intensity variance, 

equations (12) and (13) become 



(22a) a £I 2 = 1.24 k 7/6 Z 11/6 C N 2 , (Plane Wave) 

(22b) a £I 2 =0.52 k 7/6 Z 11/6 C N 2 , (Spherical Wave). 



Empirical results [Refs. 14, 32-34] have shown that when the 
optically measured o^, is plotted against Tatarski's 

theoretical values of equation (22), calculated 

from the temperature structure constant, the measured 
scintillations reach their maximum value for of the 

order of 1 to 3, then decrease slowly and tend to level off 
with increasing This approximately constant value 

for > 1 is called the saturation effect. Yura 

[Ref. 17] states that the log-amplitude covariance function 
differs both quantitatively and qualitatively in the 
saturation region from that under unsaturated conditions. 
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The actual value of at which saturation begins 

has not yet been strictly established. Strohbehn [Ref. 2] 

indicates that saturation begins for > 0.8, while 

de Wolf [Ref. 18] and Lutomirski and Yura [Ref. 1] maintain 

2 

it occurs for > 1. 

2 . Theory 

De Wolf has defined the saturation region in terms 

, 2 
of o n as 



(23) 



? k V 

1 < a £I “ — 



or equivalently, 



(24) 



1 < K m 1/2 kZ 2 C n 2 « 



De Wolf [Ref. 18] improves upon the Rytov approximation by 
introducing a directional correction factor, cos [K* p_(Z) ] , 
to allow rays being refracted by a turbule to be deflected 
by more than one Fresnel zone size off the optical axis. 

This directional factor eliminates the focusing effect of 
turbules close to the detection point and introduces serious 
interference effects and therefore amplitude changes. 

De Wolf has calculated the log-amplitude variance in 
the saturation region based on plane wave theory to be 



(25) 



< v 2 _ CK 7/3 3 2 -1/6 „ = 2tt 

Xi > " °- 41 ( K m Z C n } ’ K m > 



n 



m 
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and derives the result that the irradiance distribution 
remains log-normal in the saturation regime. However, to 
verify this equation, K m must be determined empirically, 
which is difficult to do. 

Other authors have attempted to describe the 
saturation effect by carrying higher order terms through the 
wave equation. The efforts have been extremely complex 
mathematically and have resulted in little success. 

D. TEMPORAL- FREQUENCY SPECTRUM ANALYSIS 
1. Description 

Like any other time-varying signal, a fluctuating 
voltage response from the output of a photodetector possesses 
an autocorrelation function. The power spectrum of the 
signal can be found (1) manually by computing the Fourier 
transform of the signal's autocorrelation function, or 
(2) electronically by feeding the fluctuating voltage signal 
into a Spectrum Analyzer, which performs the transform. The 
power spectrum shows the frequencies composing the signal 
and their relative magnitudes. Current theories (described 
below) indicate that the maximum frequency that will be 
present will be directly proportional to the mean wind speed 
normal to the propagation path and inversely proportional 
to the square root of the wavelength used and distance of 
propagation. Theories differ significantly on the effects 
of aperture diameter on the power spectrum. 
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2. Theories 



a. Lutomirski and Yura 

It is assumed that all refractive index 
inhomogeneities of the turbulent medium move convectively 
with mean wind speed v everywhere along the propagation 
path. Intensity fluctuations will be dominated by the 
velocity component normal to the propagation path, v . If 
the transverse intensity correlation length is p , then the 
correlation period of the intensity fluctuations at a point 
is of the order of ~ p„/v , giving a characteristic 
frequency of 



(26) 



v 

£ = EL 

o 2 tt p . 



1/2 

Using Lutomirski and Yura's [Ref. 1] model of p^ ~ (Z/k) ' , 
the maximum frequency cutoff for scintillation effects at a 
point is given by 



(27) 



v 

n 



(2ttAZ) 1/2 



The normalized photodetector signal at time t 



is given by 



C28) S’(t) ■ S(t) <S(t ' )> tj - > - 3 



ICPj.t) - <I> j2 

<T> 



d P 1 » 



and from the frozen-in turbulence assumption. 
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(29) S’ (t + x) = 



/ 



1 (P.2 ‘ v n T “ Kl> 



<T> 




-2 * 



The autocorrelation function of S'(t) is 



(30) R ,(t) = <s'(t)s'(t+T)> = -L- 

<I> 



/ 




The power spectrum is the Fourier transform of R s ,(x) and 
is given by 



Aperture averaging effects on the spectrum are seen through 



the Wq term, = -q— , where D is the aperture diameter. 

is a first-order Bessel function and C^(K) is the two- 
dimensional power spectrum of the intensity fluctuations. 
The square bracketed term falls off rapidly for w > 
displaying the filtering action of the aperture diameter. 



$ n (f), versus log f, a change in X, Z, v^ or D effects a 
shift of $ n (f) along the log f axis, 
b. Clifford [Ref. 19] 

A theoretical solution for the temporal- 
frequency spectrum of scintillation at a vanishingly small 
aperture is given as shown for plane and spherical waves 
in Figure 3. The corresponding asymptotic forms are given 



(31) $ s ,(u0 = 




2v 



When plotting a normalized power spectrum, 
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for plane waves as 



(32) <Hf) = 0.44 (o ZA 2 ) (^-) (1 + 0.27 ft 4/3 ) , 

* o 

(33) $(f) = 1.44 (c^/Hf-) «" 8/3 , ft » 1 

o 

and for spherical waves as 

(34) $(f) = 0.245 Ca £A 2 ) (f~) [1 + 0.119 ft 4/3 ] 

o 

(35) $(f) = 2.81 (° £A 2 )(f-) fl ’ 8/3 , « » 1 

o 

where ft = and f = v (2 itAZ) ^ ^ . 

±„ o n 



ft << 1 , 



9 



, ft << 1 , 



9 
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II. NATURE OF THE PROBLEM 



A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 
1 . Description 

A photodetector placed at the image point of a 
lens illuminated with a plane wave responds to the square 
of the integrated amplitude over the convergent wavefronts. 
Having passed through a turbulent medium, the wavefront of 
incoming light for an ideal spherical or plane wave is no 
longer coherent nor of equal amplitude across the surface 
at the detector aperture. The random-sized turbules 
of different refractive indices along the propagation path 
have perturbed the wavefront creating a buildup of amplitude 
in some areas and a cancellation in others. A transverse 
correlation distance, p^, is defined as the distance in which 
intensity fluctuations remain correlated on the traveling 
wavefront at the point of observation, and in general is 
a function of path length Z and wavelength of the source. 

If a detector of diameter D is used to measure the intensity 

fluctuations at a point distant from the source, then the 

* 

variance of the log-amplitude measured will depend on the 
value of D as compared to p^. For D > p^, the wavefront 
incident upon the detector will not be strongly correlated, 
cancellations from the amplitudes of different parts of the 
wavefront being out of phase will occur, and intensity 
fluctuations measured will be less than for a correlated 
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wavefront, i.e., D <_ p^. This observed decrease in the 
variance of the log-intensity fluctuations due to increasing 
aperture size is called the aperture averaging effect. 

2 . Theories 

a. Lutomirski and Yura [Ref. 1] 

For small scintillations and homogeneous 
turbulence, the transverse correlation distance of the 
intensity fluctuations is given as 

(36) p £ ~ (Z/k) 1/2 = .399 (AZ) 1/2 . 

For strong scintillation, no reliable estimate of p^ exists. 
If the aperture diameter D < p^, then the intensities 
measured at all points r^ within the aperture, I(r^,t), are 
correlated at time t, and the photodetector current s(t) 
is proportional to 

ttD 2 

(37) s(t) - 1(0, t) , 



i.e., an effective point detector. As D becomes greater 
than p^, s(t) fluctuates less and for D >> p^, the average 
incident intensity becomes constant and we have 



(38) 



<s> = 



ttD‘ 

T 



<I> 



The signal variance is given by 
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(39) 





C I I — 1 ~—2 I ^ d 




2 



where Cj-(p) is the covariance of irradiance. Lutomirski 

and Yura state that aperture averaging "...is the ratio of 
2 

a s in (equation 39) to what the variance would have been 
if the intensity over the whole aperture varied as it did 
at a point." At a point, 



(40) 




r, 2 2 

C^) 



Cj CO) . 



The aperture averaging factor. A, is defined as: 



(41) 





C l(|Pl- p 2 D 2 ,2 

— cYun — d E i d p 2 • 



Reference 1 shows that if: 

>> D, then A = 1, i.e., there is no aperture averaging 

factor ; 

P* « D, then A ~ 4 (Ji) 2 ; 

p^ ~ D, then the value of A requires a detailed knowledge 
of Cj(p). Cj(p) varies with source diameter, propagation 
slant path and meteorological conditions. For homogeneous 
turbulence and plane wave propagation, Cj(p) is not known 
for strong scintillation. It is recommended that 



31 



(42) 



A 



1 



1 + 




be used for homogeneous turbulence when the scintillation is 
small . 

b. Fried [Ref. 21] 

Using the statistics of propagation of an 
infinite plane wave over a horizontal path, Fried derives 
equation (37) , then develops an expression for the signal 
variance , 

D 

(43) o s 2 = 2 tt J pK o (p,D) Cj(p)dp , 

o 



where K Q (p,D) is the area of overlap of two circles of 
diameter D whose centers are displaced a distance p, and 
is given by 



(44) K 0 (p,D) = 



D r “ 1 r P n 

T tcos V 

0 , 



- c£M 1/2 ) , p £ D 






p > D. 



Fried then derives the same aperture averaging factor as 

given in equation (41) but labels the factor 0. He states 

that the normalized log-amplitude covariance, C^(p)/C^(0), 

- 1/2 

is a function of p(4Z/k) only, for the horizontal 

propagation case, and gives for an aperture diameter 
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normalization factor, 

(45) (4Z/k) 1/2 = .798 (AZ) 1/2 . 

Fried's theoretically calculated dependence of 

1/2 

0 upon the normalized aperture diameter, D/ (4Z/k) , is 

shown in Figure 4 for various values of 0^(0). 

c. Tatarski [Ref. 20] 

This theory on aperture averaging effects 

gives the amplitude correlation distance as being "on the 

1/2 

order of" the first Fresnel zone size, ( X Z ) . For this 

model then, the aperture normalization factor is simply 
1/2 

(AZ) . Tatarski was the first to define the aperture 

_ 2 

0 (fM Q 

averaging factor G(R) = qjQy > where Q(R) = ^ anc * 

a 2 
s 

Q(0) = — f which Lutomirski and Yura have labeled "A" 

<s> z 

and Fried "0." 

3 . Published Experimental Data 

a. Homstad and Heneghan [Ref. 22] 

Measurements were made on a mesa near Boulder, 

Colorado, over a horizontal (statistically homogeneous) path 

at distances of 200-1600m using a He-Ne (6328 A) laser. 

1/2 

G(R) vs. R/(AZ) was compared with Fried's theoretical 
curve of Figure 4. G(R) is the aperture averaging factor, 
equivalent to A in equation (41) . The result was that good 
agreement with theoretical calculations was observed 
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at shorter 



"...in the region of weak fluctuations (i.e., 
path lengths) ... data in the region of strong fluctuations 

1/7 

...show a trend toward larger G(R) for R/(XZ) > 1.0, in 
agreement with Kerr" [Refs. 23 and 24]. This is consistent 
with recent theoretical predictions [Refs. 25 and 26] and 
with Kerr's experimental work which found the correlation 
radius to be longer in the region of strong fluctuations 
[Refs. 23, 24, 35 and 36], 

Homstad and Heneghan point out that a reason 
that experimental data collected by Fried [Ref. 16] did not 
support his theoretical predictions may have been that the 
apertures he used for the 8 km path length were not large 
enough to extend beyond the correlation radius and hence 
gave a relatively constant G(R) as reported by Fried. 

b. Strohbehn and Homstad [Ref. 27] 

In an experiment conducted over land, good 
agreement with Fried's model in the region of weak fluctua- 
tions is confirmed. 

c. Kerr [Ref. 23] 

States that the predicted value of the aperture 

averaging factor should vary as D X for diameters much 

1 / 2 

larger than (AZ) , where different authors have assigned 
values of 2, 7/3 or 3 to x, and that he did not observe this. 

d. Previous NPS Work [Ref. 28] 

Schroeder observed that "Significant aperture 
averaging was observed in an enclosed laboratory for 
apertures less than the first Fresnel zone size." From 
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two field experiments (both conducted at the same path 
length and with the same wavelength) Schroeder stated that 
the aperture averaging factor may be dependent on C N and 
that "...a 4.8 mm aperture can approximate a vanishingly 
small receiver and should be used whenever possible." The 
two field experiments were conducted transmitting across 
Monterey Bay with receiver and source located on shore. 
Path length was 4.05 km, 98 percent of which was over the 
ocean surface at 4.9 meters. Aperture sizes used varied 
from .3 cm to 2.3 cm. By graphing M<cr £^( D ) “ a £ A (23mm)> 
vs. Aperture Diameter D," Schroeder* s data indicated that 
the largest average difference in values was seen for 
a 4.8 mm aperture. 

B. APERTURE AVERAGING EFFECTS ON THE TEMPORAL -FREQUENCY 
POWER SPECTRUM 

1 . Tatarski's Theoretical Model 

With increase in aperture size, Tatarski predicts 
a loss in the high frequency components of the normalized 
spectrum. By the addition of the factor 

2 

2Jj (rVk Z + 4ir 2 f 2 /v n ‘ ! 

~ R VK”Vf Z /v ' . 

“• n 

in the kernal of the integral for the plane wave model of 
W(f) , the spectral density of intensity fluctuations, a 
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partial suppression of the normalized intensity fluctuations 

as seen by an aperture of radius R is accounted for. For 
v 

f >> > this correction factor is a very small number, 

which corresponds to a suppression of high frequencies. 

As the aperture size is increased, the frequency at which 
the suppression occurs decreases. The low frequency portion 
of the spectrum will remain unaffected [Ref. 20]. 

Using the frozen-in turbulence model, the relation- 
ship between large spatial frequencies, K, and large temporal 
frequencies is given by to = K*v n . This implies that a 

receiver of diameter D will pass only temporal frequencies 




2 . Published Data 

a. Hohn [Ref. 29] 

A He-Ne (632 8 A) source was used over path lengths 

of 4.5 km and 14.5 km. Various aperture sizes from .5 cm 

to 8 cm were used. Wind speed values were obtained from 

a meteorological center near the test site. Normalized 

spectrum curves were deduced from the data taken and the 

results were that "No systematic correlations of U(f)... 

or f^/2 Were found with changing aperture diameter D of the 

photometer." A second result was that the theoretical 

- 1/2 

relation f “ Z did not appear to hold empirically. 

b. Schroeder's NPS Work [Ref. 28] 

An aperture averaging frequency analysis sequence 
was conducted during each of the two field experiments 
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previously mentioned. Schroeder noted that changes in wind 
velocity could produce a different spectrum if two data 
samples at different aperture sizes were not conducted 
simultaneously. General wind speed conditions were measured 
at the transmitting location and assumed to be the same 
over the propagation path. The spectrum was obtained by 
passing a recorded scintillation signal through a wave 
analyzer with a bandwidth of 10 Hz at 10 Hz intervals from 
20 to 100 Hz, then at 100 Hz intervals from 100 to 1000 Hz. 
The power was determined by integrating the output curves 
with a planimeter. Results were that for f < 90 Hz, the 
aperture averaging effect on the relative power spectrum 
was ’’unpredictable , " while for f > 90 Hz, the smaller 
apertures show a larger contribution to relative power, 
as predicted by Tatarski's theoretical model. A second 
result is that for the same aperture size, the power spectrum 
relative to the power at 20 Hz measured on a day when wind 
speed was estimated to be 16 knots over the path, was overall 
higher in magnitude at all frequencies above 20 Hz than for 
the day with wind speed estimated at one knot. 



C. 



J N0 



vs 



J NT 



As stated previously, one of the goals of the project 
is to be able to predict, through meteorological measure- 
ments, optical propagation conditions. A common indicator 
of this is the refractive index structure constant, , 
which is related to the log-amplitude variance through 
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equations (12) and (13). can be computed from the 

temperature structure constant, C^, [equation (6)], through 
the relationship 

(46) C T 2 = (79 x 10' 6 3p- ) -2 C NT 2 , 

where is the refractive index structure constant pre- 

dicted from temperature fluctuations. Agreement between 
the structure constant C^q obtained from long-path measure- 
ments of the log- intensity variance with may be taken 
as an indication that the optical fluctuations are not 
saturated and that the point-measured temperature fluctua- 
tions are characteristic of the path. 
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III. EXPERIMENTAL PROCEDURE 



A. GENERAL EQUIPMENT DESCRIPTION 
1 . Detectors 

To compensate for the time dependence of the 

variance, in measurement of effects of aperture size, a 

United Detector Technology PIN-5DP photodetector (Noise 

-13 

Equivalent Power 5 x 10 watts) with fixed aperture size 
was used as a reference. This was operated without bias 
voltage with the circuit shown in Figure 5. 

For the variable aperture detector, a TIXL69 
silicon avalanche photodiode (operated at its optimum S/N 
reverse bias voltage of 159V) was used to provide the high 
sensitivity at low S/N needed for measurement at very small 
aperture sizes. The characteristics of the TIXL69 are 
given in Appendix B. 

For the first series of experiments conducted in the 
enclosed corridor laboratory, both detectors operated without 
field lens systems. The field lens images the entrance 
aperture on the detector so that the light entering the 
aperture is spread over the surface of the detector. This 
avoids fluctuations due to the image wandering over the 
detector active area which may vary in sensitivity [Ref. 28]. 
For the field experiments, the following field lens systems 
were installed: 
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(a) Avalanche Detector (Figure 6) . A 50 cm focal 
length objective lens with a diameter of 4.9 cm focused the 
beam onto a 1.5 cm focal length, 9.5 mm microscope objective 
lens. Between the lenses a clear glass beam splitter 
reflected a portion of the image into an eyepiece for align- 
ment purposes. The field lens magnification was designed 

so that the image of the aperture filled the detector active 
area (1.8 x 10 cm ) . A 4.9 cm, 30 A filter centered at 

O 

6328 A with a 77 percent peak transmission coefficient 
preceded the lens system, and an 18 leaf iris, adjustable 
from 1.5 mm to 42 mm aperture diameter, was placed in front 
of the filter. 

(b) PIN-5DP Photodiode Detector (Figure 7). This 

field lens system is fundamentally the same as that for the 

avalanche detector with the exception that a 33.3 cm focal 

length, 2.54 cm diameter objective lens and 1.8 cm focal 

length, 9 mm diameter field lens were used to cover the 

2 

detector active area of .05 cm . 

2 . Fan-Scan Mode of Laser Operation 

O 

A CW Radiation 5 mWatt 6328 A laser operating in the 
TEM oo mode was used for the aperture averaging experiments. 
The NPS "Fan-Scan" transmitter optics developed for ship- 
to-shore measurements was used. This system broadens the 
circular beam cross-section horizontally by reflection from 
a variable-radius cylindrical mirror, then scans this "fan" 
vertically by reflection from an oscillating output mirror. 
Thus the narrow uniform beam scans vertically twice each 
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cycle over a detector within a square field of view of order 
20 meters side at 1500 meters (variable). Ship motion and 
background compensation are provided by measuring the dif- 
ference between maximum and minimum irradiance values. 

3. Electronic Components 

Figure 8 is a block diagram of the electronic 
components used in all aperture averaging and frequency 
analysis experiments. The voltage signal from each detector 
package was sent to a PAR 113 (Princeton Applied Research) 
low noise differential pre-amplifier with maximum gain of 
10 K. The outputs from the two PAR'S were connected to a 
single switch box, from which the selected output passed 
to the demodulator. This switch box enabled the operator 
to rapidly record alternate scintillation signals from the 
two detectors. The demodulator stretched the voltage peak 
and sampled the stretched value, sampled the noise voltage 
100 microseconds later in between peaks, and sent the dif- 
ference to the Log Converter. The log-converted difference 
was then sorted by the PIP into the appropriate voltage 
channel as one count. The demodulator provides a trigger 
pulse for a signal monitoring oscilloscope and the PIP. 

The calibration box provides calibration voltages 
of .5, 1, 2, 4 and 8 volts to the Log Voltmeter/Converter, 
and from there into the appropriate voltage channel on the 
PIP 400 Multichannel Pulse Height Analyzer (MPHA) . A 
readout of the counts in each channel during calibration 
is used by the HP9810A Calculator and Plotter to scale the 
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voltage. The PIP 400 MPHA records the number of times a 
particular voltage value (i.e., log- intensity) is measured 
from the fluctuating signal during a run. This information 
is displayed durin'g collection and can be stored indefinitely 
for subsequent printout and analysis as long as no other 
data is recorded further on that channel. The scintillation 
signal for frequency analysis was recorded from the output 
of the Log Voltmeter/Converter . 

4. Shipboard Detector Stabilization 

The avalanche and PIN photodiode detectors were 
mounted on a feedback-controlled gyros tabilized platform. 

The system is designed to compensate for ship motion and 
keep the critical detector alignment stable. Currently this 
system’s performance is adequate at short ranges (less than 
2000 m) and fairly calm sea states, but is degraded rapidly 
by increased distances or moderate sea states. In these 
instances the detector (gyroscope) alignment must be con- 
trolled manually by keeping the laser signal image in the 
center of two crosshairs of an eyepiece mounted on the 
gyro support. Gyroscope nutation, inability to keep the 
detectors aligned properly, and movements of the ship all 
contribute errors to a measured value of signal variance. 

B. APERTURE AVERAGING EXPERIMENTS 

1 . Propagation Paths and Turbulence Conditions Present 
a. Corridor Laboratory 

Aperture averaging experiments were conducted 
in a long basement corridor by locating a 5 milliwatt He-Ne 
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(6328 A) laser at one end of the corridor, propagating the 
beam the length of it, reflecting it back down the length 
of the corridor using a plane mirror and measuring the 
scintillation near' the location of the source. Total path 
length was 280 m. The avalanche and PIN photodiode 
detectors were clamped to a vertical support. Once aligned, 
the detectors were not subject to any of the inherent motions 
from platform instability which are present when conducting 
the same experiment on board a ship. The experiments were 
conducted during early morning hours (with the exception 
of that on 7 April) to eliminate rapid variances in turbu- 
lence conditions due to persons opening and closing exterior 
doors. Turbulence was provided by ventilation system openings 
and building heating ducts located near the propagation path. 
Earlier measurements on the nature of the turbulence in 
this corridor indicated that the turbulence present did not 
represent homogeneous (Kolmogorov) turbulence conditions. 
However, the non-homogeneous turbulence did produce an 
intensity distribution that proved to be log-normal as 
predicted from theory. A sample log-intensity distribution 
taken during a corridor experiment is shown in Figure -9. 
b. Over Land 

To complete an aperture averaging sequence in 
homogeneous (Kolmogorov) turbulence before going to sea, 
where the ship’s motion can affect the measured signal 
variance, an aperture averaging sequence was conducted on 
the Naval Postgraduate School Annex grounds on which the 
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Fleet Numerical Weather Central is located. The site 
selected to conduct the experiment consisted of a level 
terrain, three fourths of which was a blacktop covered road, 
the remainder a level dirt road. There were buildings and 
trees within 50 m of the propagation path; however, no 
forced draft blowers, heating systems or machinery was 
present to disturb the apparent homogeneity of the path. 

Winds were negligible; hence no eddy effects from winds 
curling around buildings or trees were present. The 5 mW 

O 

He-Ne (6328 A) source was mounted on a large telescope 
frame in a small, open-end trailer. It transmitted a 
Fan-Scan signal over a path length of 336 m to the avalanche 
and PIN photodiode detectors mounted on a metal plate 
affixed to a manually trained telescope. The equipment 
setup of Figure 8 was used and a frequency analysis experi- 
ment was conducted in conjunction with the aperture averaging 
Again, intensity distributions over the path proved to be 
log-normal . 

c. Over Water (Land to Land) 

Two experiments were conducted over different 
path lengths and under different turbulence conditions. 

Experiment I: Coast Guard Wharf to Spanagel 

Roof. The optical path, of length 2335 m, was highly 
irregular. The propagation path was slanted downwards, 
going from the roof at an altitude of about 75 m above sea 
level down to the Mobile Optical Research Laboratory at the 
Coast Guard Station parking lot at about 10 m above sea level 
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About one third of the path was over water, the remainder 
over several commercial and residential structures. During 
the course of the experiment there were intermittent rain 
squalls. Data revealed that this entire sequence was 
conducted under saturation conditions for . 

Experiment II: Monterey Holiday Inn to Hopkins 

Marine Station. The second experiment, conducted over a 
propagation path length of 4230 m, was 97 percent through 
the marine boundary layer. The laser beam was transmitted 
over a horizontal path at about 10 m above sea level. The 
sky was clear on 11 March with average wind speeds during 
the experiment of about 7 m/sec. The He-Ne (6328 A) laser 
source was located in the rear of the open-end trailer and 
was operated in the Fan-Scan mode. The trailer was located 
at the Monterey Holiday Inn parking lot, about 100 m from 
the water's edge. The terrain between the ocean and the 
trailer consisted of a low, sandy incline sloping downwards 
from the trailer's position towards the ocean. The detection 
devices and processing system were located in the bus (Mobile 
Optical Research Laboratory) stationed on the Hopkins Marine 
Station grounds about 10 m above the ocean surface and 20 m 
from the water's edge with an unobstructed view of the 
signal source. 

Experiments measuring MTF being conducted 
simultaneously with aperture averaging precluded the PIN 
photodiode detector's being available for use as the fixed 
aperture system on this occasion. For this experiment a 
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HgCdTe detector with an aperture size of 23.6 mm was used 

as the fixed aperture system. The HgCdTe detector measured 

scintillation from a (10.6 pm) CO^ laser. It is assumed 
2 

that measured at 10.6 pm varied temporally in the same 

way as measured at .6328 pm, as indicated by the consistency 

2 

of multiwavelength measurements of . 
d. Over Water (Land to Ship) 

The He-Ne Fan- Scan transmitter on the bus was 
located at Pt. Pinos (Figure 10) as near to the shoreline 
as possible (about 50 m away) , with the detectors and data 
processing system on board the Research Vessel ACANIA. Only 
those runs were recorded in which detector alignment was 
maintained for periods of at least, two to three minutes. 

Path lengths used during these runs were 1442 m, 1550 m, 
and 1625 m; the maximum range was limited by the failure 
of the tracking system to lock-on automatically at ranges 
greater than 2000 m. Path lengths were determined by base- 
line triangulation methods on shore using surveyor's 
transits. The optical paths were about 98 percent over the 
ocean surface. 

2 . Data Reduction 

Data reduction, i.e., determining the refractive 
index structure constant and the log- intens ity variance 
from a fluctuating voltage signal, was accomplished within 
minutes of the completion of a data collection run through 
the use of the Victoreen PIP 400 Multichannel Pulse Height 
Analyzer and the Hewlett Packard 9810A Calculator and 
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Plotter. As described previously, the PIP 400 stored the 
log-intensity distribution of the irradiance measured. 

An interface was developed between the PIP 400 and the 
HP9810A to allow the calculator/plotter to print out the 
log- intensity distribution and determine and o ^ from 
it. The algorithm the computer followed [Ref. 31] was to 
apply a best- fit Gaussian curve to the distribution to 
obtain , and then determine through the relation 

0 o i 

C4?) ^ = k" 1 * Z ii/1,r 'c.4763 i/2 ' 

An example of the printout is shown in Figure 9. The three 
markers near channels 60, 116 and 173 are calibration 
points of .5 V, 2 V and 8 V respectively, which are equally 
spaced on the log scale. 

3. Aperture Averaging Factor G(D) 
a. Description 

To determine the aperture size that best 
approximates a point receiver, i.e., that size at which 
reaches a maximum for a given path length and 
value and does not decrease significantly for further 
reduction of aperture size, the following aperture averaging 
factor was used: 
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(48) G(D) 



C N f V ( D o > VI C N t F (D o ’ C N [V ^ D o ’ V 1 

= c^L'FC^/tojJ * Cjfivnr,tjT ' “ q7lT0)*,tTJ" • 

Cj^ [V(D Q , t Q ) ] is the value of measured using the variable 
aperture detector at the maximum diameter D Q and the 
time t . Similarly, F denotes the fixed aperture detector 
which is used to measure the change in due to variables 
other than a change in aperture size. The first factor 
C N [V(D 0 ,t 0 )]/C N [F(D 0 ,t 0 ) ] , compares the initial C N value 
measured (at time t„) with the variable diameter set at D 
to that measured by the fixed diameter at time t . 

[F (D Q , t) ] /C^ [V(D, t) ] reflects the change in that ratio 
produced by a reduction in aperture diameter from D q to D. 

If the two apertures (fixed and variable) measure simul- 
taneously at times t and t, and the variable aperture size 
is not changed between the two measurements, then G(D) 
will ideally equal 1.00. This was tested in conjunction 
with two corridor experiments conducted on 10 February 
and 5 March 1976. The comparison runs produced G(D) values 
of .9818 and .9981, indicating an average error of about 
one percent for G(D) values at different aperture sizes 
due to experimental variances. 

b. Reason for Use; Relation to Fried's 0(D) 

By starting with the largest experimental 
variable aperture diameter, D Q , for which G(D q ) = 1.00, 
then progressively reducing it until the minimum variable 
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aperture D m ^ n is reached, the aperture size that best approxi- 
mates a point detector can be determined by noting at which 
aperture diameter the lowest value of G(D) occurs. Empirical 
results (Section IV. A.) show that G(D) does not remain linear 
from D down to the smallest possible aperture size, D . , 

but levels off and reaches a minimum value, [G(D)] ■ , well 

’ L J min 

before it. 

Fried [Ref. 21] defined an aperture averaging 
factor 0(D) [equivalent to A in equation (41)] based upon 
a vanishingly small receiver (D -* 0 mm). Fried's theory 
states that all apertures of diameter D > 0 will measure a 
signal variance of smaller magnitude than would be measured 
using a vanishingly small detector. When using an aperture 
of diameter D, the measured variance must be multiplied by 
the correction factor to allow for aperture averaging. 

[G(D) ] • identifies the aperture diameter, D . , that best 
approximates a vanishingly small receiver. G(D) is related 
to Fried's aperture averaging factor in the following manner: 

G(D) StVCD.Ol C N [V(D 0> t)] _ x 

(49) Ttrr Ttrr 

C. APERTURE AVERAGING EFFECTS ON TEMPORAL -FREQUENCY 
SPECTRUM ANALYSIS 

1 . Experimental Conditions 

a. (7-9) April 1976, Basement Corridor 

The electronic equipment setup used in aperture 

averaging sequences is shown in Figure 8 and the procedure 
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followed in obtaining and processing signals on this 
date is described in Sections III.B.l.a, IV. A. 1, and 
IV. A. 2. a. (3) . For frequency analysis, the scintillation 
signal was recorded from the output of the Log Voltmeter/ 
Converter with a Precision Instruments, Model No. PI-6204, 
four-channel recorder at 37.5 ips. Performance characteris- 
tics of the recorder are contained in Appendix C. Signals 
were recorded for a period of about 75 seconds each at 
aperture sizes of 40, 30, 20, 15, 10, 7, 5, 3, 2 and 1.5 
millimeters . 

b. 20 April 1976, Navy Annex 

Experimental conditions are as described in 
Sections III.B.l.b, IV. A. 1, and IV.A.2.C. Frequency analysis 
data was recorded as described above for the (7-9) April run. 
Aperture diameters examined were of 40, 30, 20, 15, 10, 7, 

5 and 3 millimeters. Power to operate the equipment on 
the bus was provided by a portable gasoline-powered motor 
generator. 

c. (29-30) April 1976, R/V ACANIA 

Experimental conditions on this date and 

procedures followed for aperture averaging are described 
in Sections III.B.l.d, IV. A. 1, and IV.A.3.(d and e) . The 
signal to be frequency analyzed was recorded from the 
output of the Log Voltmeter/Converter (Figure 8) on a 
Sanborn HP- 3914, 14-channel Analog Tape Recorder, operated 
by the meteorology team recording micrometeorological data 
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for comparison. Operating characteristics of the 
HP-3914 Recorder are included in Appendix C. Time reference 
signals were recorded on the tape for identification. 

2 . Data Reduction 

The recorded scintillation signal was frequency 
analyzed by playing back the signal into an EMR Schlumberger 
Instruments 1510 Digital Spectrum Analyzer. The spectrum 
analyzer was operated in the continuous mode over a frequency 
range of 10 Hz to 2.56 KHz. The signal was sampled by the 
spectrum analyzer for a period of 25.6 seconds. The resulting 
spectrum was plotted on a Hewlett Packard Moseley 7035 B X-Y 
Recorder. Frequency was plotted linearly along the X-axis 
from 42 Hz to 2.56 KHz. The Y-axis was a 0 to -60 dB linear 
scale . 

For evaluation of the aperture averaging effects on 
the spectral density, relative power spectra were obtained 
directly by aligning the -60 dB line of the spectrum of each 
different aperture to a common new reference axis and 
plotting each spectrum manually. Points of local minima 
were plotted on each spectrum, thereby avoiding abnormal 
magnitudes due to noise spikes. 
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IV. EXPERIMENTAL RESULTS 



A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 

1 . Commonalities 

In addition to the electronic equipment setup 
(Figure 8), detectors, and laser source previously described, 
the following items were common to all aperture averaging 
experiments: (1) The Hewlett Packard 9810A Calculator and 

Plotter was used for immediate data reduction. (2) Aperture 
sizes for the variable aperture system were set by inserting 
a micrometer caliper set at the proper diameter size within 
the aperture and closing the iris upon it. Both horizontal 
and vertical dimensions were checked. The aperture size was 
rechecked at the end of the run. 

2 . Preliminary Investigations 
a. Corridor Experiments 

(1) 2235, 7 February 1976 - 0515, 8 February 1976 
This initial aperture averaging sequence was 
conducted to examine the validity of the assertions made 
by Schroeder [Ref. 28] that a 4.8 mm aperture diameter best 
approximates a vanishingly small receiver. For this experi- 
ment, the laser was operated in the Scan mode only and a 
neutral density filter of .5 was used to decrease the signal 
strength. The field lens detector systems of Figures 6 and 
7 had not been constructed at this time. Each detector 
consisted of the photodiode placed at the focal plane of an 
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objective lens in a narrow cylinder of diameter 22.6 mm with 
a 10 A filter centered at 6328 A at the cylinder front. This 
made the largest variable aperture size available 22.6 mm. 
Aperture sizes of '22.6, 20, 17, 14, 12, 10, 7, 6, 5, 4, 3 
and 2 millimeters were used. G(D) values calculated from 
the measured log-intensity distribution standard deviations, 
a^j, are given in Table I with the plotted "G(D) vs. Aperture 
Diameter, D M shown in Figure 11. The last column of Table I 



gives 

averaging. 



G(D) 

TO! 



mm 



the correction factor for aperture 



In Figure 11, the curve bottoms out at an 
aperture size of about 5 mm and the corresponding value of 
G(D) does not significantly decrease from its value at 
that point as the aperture diameter is further reduced. 

(2) 1655-1820 and 2100-2240, 10 February 1976 
It was noted from the experiment of 
7 February that although the log-intensity distributions 
were for the most part Gaussian, on several of the larger 
aperture runs double peaks occurred near the top of the 
log-intensity curve. It was believed that this was due to 
the laser being in the Scan mode only so that the detector 
would see one average intensity peak as the beam swept up 
and a slightly different one as the beam swept down. In an 
attempt to alleviate this, a second aperture averaging 
sequence was conducted with the laser in the Fan-Scan mode. 
By spreading the beam out and scanning it, it was felt that 
a more uniform intensity might be seen by the detector, and 
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the double peaks eliminated. A Fan-Scan divergence setting 

of seven milliradians was used and a neutral density filter 

of optical density .5 was inserted in the beam to reduce the 

signal intensity at the detector. With the beam reflected 
2 

off a 625 cm mirror at a distance of 140 m, the estimated 

- 9 

signal strength at the position of the detector is 1.7 x 10 
Watts/cm , which is of the order of the signal intensity 
for the at-sea experiments, described later, at a path 
length of 1600 m. 

Aperture sizes used were identical to those 
used during the 7 February run. Log Intensity distributions 
on this day revealed that the double peaks were eliminated 
in all but one distribution by use of the Fan-Scan, and 
this mode was therefore adopted as the standard laser opera- 
tion mode for future aperture averaging sequences. Figure 
12 shows "G(D) vs. Aperture Diameter, D M for this experiment. 
Although this experiment was conducted identically to that 
on 7 February with the exception of the Fan-Scan operation, 
three significant differences stand out in the results. 

First, the curve bottoms out at the 6 mm aperture size vice 
the 5 mm size of 7 February. Secondly, G(D) values at all 
aperture sizes are in general larger than they were on 
7 February, even though the same value of D q ( 22.6 mm) was 
used to calculate G(D) values. Thirdly, the shape of the 
curve following the bottoming out at 6 mm is different from 
that seen on 7 February. In comparing the data, the 
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average [F(D q , t) ] is seen to be significantly different 



on the two days, implying a larger C^. On 7 February 
< ^l[F(D 0 .t)] > = .330, while on 10 February <a ^ ^ [F (D , t) ] > =.428 
An increasing <a ^j > seems to increase G(D) values, as can be 
seen by comparing Table II with Table I. 



(3) (7-9) April 1976 

An aperture averaging frequency analysis 
experiment was conducted. The aperture averaging results 
are included for comparison with runs conducted on 7 February 

O 

and 10 February. The same He-Ne (6328 A) laser was 
operated in the Fan-Scan mode with a divergence setting of 
seven milliradians and a neutral density filter of .5 was 
inserted in the beam. 

Prior to this run the largest aperture size 
available on the variable aperture detector was 22.6 mm. 

To increase the largest aperture size available and improve 
signal detection capability, the field lens systems of 
Figures 6 and 7 had been constructed and were in use for 
this experiment. The significance here is that with a 
larger D q available, one would expect to see a smaller 
[G(D)] m i n f° r this sequence. This proved accurate, as 
shown in Figure 13. A much lower [G(D)] m ^ resulted using 
the 40 mm aperture as D q vice 22.6 mm. This curve appeared 
to follow the same basic pattern as those for 7 February 
and 10 February with the exception that this took on a 
different shape from either of the previous two for apertures 
less than the bottoming out diameter of about 5 mm. In 
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TABLE II 
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.716 .426 .594 



order to see better how this run compared with the previous 

two, G(D) values were recomputed using the 20 mm diameter 

as D . The result is shown in Figure 14. fG(D)l . for 

(7-9) April is still well below that of 7 and 10 February. 

The average for the fixed aperture system for this run 

was <cr £i , t) ] > = .08 which is much lower than the average 

values for 7 and 10 February as is [G(D)] . . Table III 

mm 

gives a comparison of the average signal variance with 
[G(D)] m ^ n for the three experiments conducted in the same 
manner, with the same equipment, but with different turbu- 
lence conditions present as indicated by <cr 0T > . 



TABLE III 



Date 


<cr £ [F(D o ,t)]> 


[G (D) ] • 

1 J J mm 




A [G (D) ] . 

1 v ' J min 


10 Feb 


.428 


.647 






7 Feb 


.330 


.468 


.098 


.179 


(7-9) Apr 


.084 


. 223 


.246 


.245 


It appears 


that [G(D)1 • increases with 
L v J J mm 


<a^j> for low values 



of , and that in general, for a given D Q , the value of 

[G(D)] m in is dependent upon the level of turbulence. The 
higher the turbulence level, the higher the [G(D)] m ^ n value. 



b. Aperture Diameter Normalization 

Examination of Figure 14 shows that the shapes 
and slopes of the plots are similar from the largest aper- 
ture size down to about 5 mm where three different forms of 
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the curve then emerge. Following previous theoretical 



models [Refs. 1, 20 and 21], it is expected that D m ^ n > the 
aperture best approximating the "vanishingly small" receiver, 
should be directly related to the intensity transverse 
correlation length, p^. This quantity is given in the 
three models as: 

(1) (AZ) 1/2 , [Tatarski, Ref. 20] 

(2) (£V /2 , [Fried, Ref. 21] 

(3) ( 1/2 f [Lutomirski and Yura, Ref. 1], 
Aperture diameters for the three corridor experiments were 
scaled according to the three models, the results of which 
are shown in Figures 15, 16 and 17, respectively. The 
three curves bottom out at approximately the value of 

Z 1/2 

D = (jr) predicted by Lutomirski and Yura. 

c. Over Land Experiment (2128-2347, 20 April 1976) 

Aperture diameters examined were 40, 30, 20, 

15, 10, 7, 5 and 3 millimeters. The Fan-Scan was set at a 
divergence of 7 milliradians and no neutral density filter 
was used to decrease the signal intensity. The laser was 
situated about 2 meters and the detectors 2.5 meters off 
the ground. The sky was clear, temperature averaged 10°C 
and winds were negligible (0-2 m/sec). 

"G (D) vs. — vrr " was plotted to locate the 

(Z/k) i/Z 

bottoming out point. Although able to test only a few 
aperture sizes, previous runs conducted outside the corridor 
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laboratory (discussed in section d below) indicated that the 
1/2 

(Z/k) normalization factor was accurate in predicting 

1/2 

the bottoming out point of the graph at D/(Z/k) =1.0. 

Figure 18 shows the graph for the experiment conducted on 

20 April. The bottoming out occurs for the aperture size 

1/2 

nearest, but not less than, D/ (Z/k) = 1.0. Note, for 

1/2 

this path length (Z/k) ~ 5.8 mm, and although the 5 mm 

1/2 

aperture size is nearer to (Z/k) than the 7 mm size, 

previous runs have shown that the behavior of G(D) below 
1/2 

D/(Z/k) = 1.0 is erratic. Thus, the 7 mm size is nearer 

to the minimum G(D) value that corresponds to the bottoming 

1/2 

out point at D = (Z/k) than is the 5 mm size. As a 

further comparison, the corridor runs conducted on 

7 February and (7-9) April at a path length of 280 m were 

plotted with this run of 20 April at the path length of 336 m, 

shown in Figure 19. It appears that the bottoming out point 

1/2 

occurs very close to D/(Z/k) = 1.0. 

[G(DJ] for this experiment, with respect to 
a 20 mm D q , was .666 and the average was 

<<J £I [F (D q » t ) 3 > = *162. This does not fit in with the results 
given in Table III; however, the propagation path and turbu- 
lence conditions were different. To verify a predictable 
increase of [ G (D)] min with < ^£j > would require repetition 
of the experiment under identical experimental conditions 
and over the same propagation path where the turbulence 
would have the same characteristics. Time considerations 
precluded this procedure. 



60 



d. Over Water (Land to Land) 

(1) 2220, 5 March - 0345, 6 March 1976 (Spanagel 
Roof to Coast Guard Parking Lot) 

Due to the irregular propagation path used 
and corresponding turbulence conditions over a long distance, 
this entire sequence was conducted under saturation con- 
ditions for o 0T . <a PT [F(D ,t)]> was 1.01 and [ G (D) ] . 

o min 

was .8808 with respect to a D q of 24.3 mm. The graph of 

M G(D) vs. Aperture Diameter, D" is shown in Figure 20. 

1/2 

For this path length, (Z/k) ' =15.3 mm. The main item 

of interest for this run is that, under saturation conditions, 
the maximum difference in G(D) values observed was only 
[AG(D) ] = .167, which included comparisons of aperture 

sizes of 24.3, 20, 16, 12, 7, 6, 5, 4, 3 and 2 millimeters. 
This indicates that aperture averaging has only a small 
effect on the measured signal variance under saturation 
conditions for aperture sizes less than 24.3 mm. 

(2) 1550-1640, 11 March 1976 (Monterey Holiday 
Inn to Hopkins Marine Laboratory) 

Due to the limited time available for 

conducting the experiment on this day, only aperture sizes 

of 41.7, 24.3 and 14.0 millimeters were examined. A plot 

of M G(D) vs. Aperture Diameter, D" is shown in Figure 21. 

1/2 

The path length was 4230 m, giving a (Z/k) of 20.6 mm. 
Figure 21 contains insufficient data to draw definite 
conclusions. The same data is plotted again in Figure 24 
in terms of the normalized diameter, and will be discussed 
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at that point. [G(D)] m ^ n and <0 £j [F (D q , t) ] > for this run 
should not be compared with other experiments in this sec- 
tion since a different source and fixed aperture detector 
system was used and is wavelength dependent for identical 

turbulence conditions. 

3. Propagation in the Marine Boundary Layer (R/V ACANIA 

Experiments.) 

a. 1840-1930, 30 March 1976 

The winds during the time of this experiment 
varied from 8 to 12 m/sec in a direction of about 10 degrees 
off the normal to the propagation path. The sky was partly 
cloudy, temperature about 7°C, sea conditions were moderate 
with average wave heights of about 1 meter, and swell heights 
of about 2 to 3 meters moving in a direction normal to the 
propagation path. The motion of the R/V ACANIA at anchor 
made maintenance of gyro lock-on difficult. This undoubtedly 
caused intensity variances to be greater than would be ob- 
served normally. The path length on this day was 1600 m. 

The detectors were located about 7 meters and the source 
about 5 meters above the ocean surface. Aperture sizes 
compared were 38, 31, 24.3, 17 and 10 millimeters. For 
this path length, (Z/k) 1 ^ 2 = 12.7 mm. The M G(D) vs. 

Aperture Diameter, D M curve is shown in Figure 22. This 
is not typical of past aperture averaging results, but 
resembles those obtained under saturated conditions. How- 
ever, <o n [F(D ,t)]> for this experiment was .149, well 
below the saturation region. It is felt that the sea 
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conditions, rapid ship movements, and inability of the gyro- 
scope lock-on system to work effectively in moderate seas 
may have obscured any detectable aperture averaging effects. 

b. 1340-1415, 2 April 1976 

Wind and sea conditions on this day were 
identical to those of 30 March. The sky being clear was 
the only difference in weather conditions. The sea state 
was still moderate and there was still considerable ship 
motion. The path length on this day was 1600 m. The gyro- 
scope lock-on system had been improved slightly by readjusting 
the existing circuitry; all other experimental conditions 
were as described for 30 March. Aperture sizes used for 
this run were 40, 30, 20 and 10 mm. Figure 23 shows 
M G(D) vs. Aperture Diameter, D" and Figure 24 compares 

i/7 

"G(D) vs. D/(Z/k) ' " for this run with the data taken over 

water with a path length of 4229 m on 11 March. Once again, 

1/2 

it appears the normalization factor (Z/k) successfully 
predicts the diameter of the aperture that best approximates 
a vanishingly small receiver. 

c. 1445-1657, 27 April 1976 

The sky was clear at 1445 but became overcast 
at about 1545. Winds were 8 to 10 m/sec blowing normal 
to the propagation path; waves also moved normal to the path 
with a height of about 1 meter and swells about 2 meters 
high. The ship rolled considerably in the swells and once 
again gyroscope lock-on was difficult. The ship was 
anchored, path length was 1650 m, and aperture diameters 
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examined were 40, 35, 30, 25, 20, 16, 13, 10, 5 and 4 

millimeters. "G(D) vs. D/(Z/k)^ 2 " for this run is shown in 

Figure 25. G(D) values seem erratic; this is attributed to 

difficulty in stabilizing the detector and tracking the 

ACANIA from shore. 

d. 1830-2024, 29 April 1976 

The sky was clear, wind speed throughout the 

experiment was 3 to 4 m/sec, and path length was reduced 

to 1442 m. Sea conditions improved slightly with wave 

heights down to about .7 m and swells of about 1.5 m. The 

wave direction was about 15 degrees off the normal to the 

propagation path as was the wind. Aperture sizes of 40, 35 

30, 25, 20, 16 and 13 millimeters were used. "G(D) vs. 

1/2 

D/(Z/k) '* " is plotted in Figure 26. The 35 mm aperture 

size was used as D Q for the G(D) calculations since it 

appeared to give more credible overall values based on 

past experiments. The G(D) value for the 40 mm aperture 

is plotted relative to the 35 mm aperture as D . With the 

exception of the 40 mm aperture [ 5 — __ = 3 . 33 ] and the 

(Z/k ) 1/2 

16 mm aperture [ — 2 = 1 . 33 ], the curve appears to 

(Z/k ) i/2 

parallel results of past runs, i.e., reaches a minimum at 

^ 1 > n ~ 1.0 and shows a nearly linear decrease of G(D) 

(Z/k ) i/2 

with decreasing aperture diameter to this point. For this 
run, [G(D)] min = .721 and < 0 n [F (D Q , t) ] > = .145. 
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e. 1044-1311, 30 April 1976 

The sea state was relatively calm with wave 

heights of about .5 m and no apparent swells. The sky was 

clear, temperature about 21°C and winds were 3 to 4 m/sec. 

The calm sea state enabled the gyroscope tracking system to 

work well in the automatic mode and results more consistent 

with past sequences were obtained on this day. The ship was 

at anchor and the propagation path length was 1625 m. 

Aperture sizes examined were 40, 16, 13, 10 and 5 milli- 

1/2 

meters and (Z/k) for this path length was 12.8 mm. 



The plot is shown 


in Figure 27. 


Bottoming out occurs near 


D/(Z/k) 1/2 = 1.0, 


l^min * ' 


766, and for this run 


<0 m [ F ( D 0 » t )] > = • 


264. A comparison of [G(D)] m ^ n and 


<a 5,I > ^° r t ^ le t ^ lree R/V ACANIA 


experiments is given in 


Table IV. 


TABLE 


IV 


Date 


[ G (D) ] . 


<a w [F(D 0 ,t)]> 


2 Apr 


. 792 


. 726 


30 Apr 


.766 


.264 


29 Apr 


. 721 


.145 



From the data of Table IV it appears again that in general 
[G C D D ] min increases with • 
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4. Comparison of Theoretical Aperture Averaging Factors 
with Experimental Results 

a. Lutomirski and Yura Model [Ref. 1] 

For aperture diameters on the order of the 

intensity correlation length, as utilized during these 

experiments, Lutomirski and Yura "recommend" an aperture 

averaging factor of 



(46) 



A = 



1 * 



Table V lists Lutomirski and Yura's theoretical value of A 

computed for our experiments and the observed value of A. 

2 2 

A, theoretically defined as a /a , is equivalent to 

s s o 

our [G(D) J m ^ / [G(D q ) ] where D q replaces D in equation (46). 

The theoretical values do not correlate well with the 
empirical values. Further, Lutomirski and Yura's theoretical 
model does not contain any turbulence dependence which appears 
to be necessary in view of the apparent dependence of 
[G(D)] m i n 0n <a £i > seen empirically. 

b. Fried's Model [Ref. 21] 

To compare empirical results with Fried's 

1/7 

theoretically derived plot of "0 vs. D/ (4Z/k) " 

(Figure 4) , whose curves are based on an aperture diameter 

1/2 

of D = .l(4Z/k) ' as representing a vanishingly small 
receiver, it was necessary to use the experimental aperture 
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30 Apr/ACANIA .290 .587 



size for each run that gave the minimum G(D) value as the 

vanishingly small aperture diameter. The effect of this is 

to shift the empirical curves to the right of the theoretical 

curve. Figures 28 and 29 show Fried's theoretical plot 

for log amplitude variance 0^(0) = 2 ^ as compared with 

experimental values for various runs. For this graph, 

aperture diameters were appropriately normalized for com- 

1/2 

parison by the (4Z/k) factor. Using equation (45] , 

0 was calculated as 



(47) 



° s hv(D.t 0 n 

^Vco-.to)] ’ 



2 

where o [V(D',t )] is the signal variance measured at 

time t and with variable aperture diameter D', which 

corresponds to the aperture diameter that best approximated 

a vanishingly small receiver, i.e., corresponded to 
2 

[G(D)] m in* G s [V(D,t 0 ]] is the signal variance that would 
have been measured at time t Q had aperture size D been used, 
and was found from the relation 



(48) 



o s 2 [V(D,t 0 )] - a s 2 [V(D, t) ] 



°s^ F(D o’V 1 > 

a s 2 [F(D 0 ,t)] 



where again, F is the fixed aperture detector. 

Empirical results appear to follow Fried's 
theoretical model generally well. In contrast with 
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Lutomirski and Yura's aperture averaging model. Fried 

acknowledges the dependence of the aperture averaging factor 

on the turbulence conditions, as represented by different 

curves for particular values of C^(0) in Figure 4. Fried's 

model indicates that as C^(0) increases, 9 increases for a 

specific aperture diameter, or the inverse of the aperture 

averaging factor decreases. This is in agreement with 

empirical results noted in Sections 2. a. (3) and 3.e. above. 

From Figure 28 it is evident that the aperture averaging 

curve for 7 February lies below that of 10 February and above 

that of (7-9) April. The 9 values are therefore in 

general larger for 10 February than for 7 February, i.e., 

2 2 

the ratio of a s /a is larger, and the inverse of the 

aperture averaging factor is less, as predicted by Fried's 
theory. 

B. APERTURE AVERAGING EFFECTS ON THE RELATIVE TEMPORAL- 
FREQUENCY POWER SPECTRUM 

1 . General Results 

The spectrum analyzer frequency limits were selected 
to examine frequencies within the range 10 Hz to 2.56 KHz. 

The demodulator in the signal processing system eliminates 
frequencies above 2 KHz, coinciding with the original 
estimate of the Fan-Scan frequency. The bandwidth is further 
reduced to a maximum frequency of 1 KHz for experiments 
conducted using the Sanborn HP- 3914 Recorder. Power 
spectra from all frequency analysis experiments show the 
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actual Fan-Scan frequency to be 1.92 KHz, clearly evident 
as a noise spike in Figures 34, 35 and 36. Previous work 
at NPS by Schroeder [Ref. 28] examined frequencies from 
20 Hz to 1 KHz. 

Results of the experiments conducted on (7-9) April, 

20 April and (29-30) April, under varying turbulence 
conditions, all indicate that aperture averaging has a 
discernable effect on the observed power spectrum. Empiri- 
cal results support the theoretical prediction that a 
corresponding reduction in high frequency power results as 
aperture diameter is increased. The reduction in power, 
however, appears to occur over a large portion of the fre- 
quency range, not solely on the high frequency end. Figures 
30, 31, 32 and 33 show the effect of aperture averaging on 
the relative power spectrum. 

2 . Turbulence-Dependent Forms of the Frequency Spectrum 

As previously remarked, the three frequency analysis 
experiments were conducted at different locations and under 
different turbulence conditions. The resulting power 
spectra for each run were similarly different. Typical 
spectra for the corridor experiment (inhomogeneous turbulence) , 
the Navy Annex run (over land in homogeneous turbulence) , and 
the R/V ACANIA data (over water in homogeneous turbulence) , 
are shown in Figures 34, 35 and 36, respectively. The 
significant item of note in Figure 34 is the uniformity of 
the spectrum below .85 KHz. The prevalence of noise spikes 
at 110 Hz and 120 Hz intervals in Figure 35 is believed due 
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to the extremely noisy, gas-powered portable generator 
used to provide power to the electronic equipment during 
this experiment. Also peculiar to this experiment were 
characteristically low magnitude power spectra. Figure 36 
shows an approximately linear decrease in power with 
increasing frequency. 

Eliminating the noise spikes in Figure 35, it 
appears that the spectrum assumes the same general form as 
that of Figure 34. A significant difference in the experi- 
mental conditions for 9 April and 20 April as opposed to 
those of 30 April was the average wind speed present. 

On 9 April and 20 April winds were negligible. On 30 April, 
conducted at sea on board the ACANIA, average wind speeds 
were about 4 m/sec. Further experiments need to be con- 
ducted to examine the effect of wind speed on the shape of 
the power spectrum. 

3. Comparison with Clifford’s Theoretical Model 

Clifford's theoretical model for the temporal- 
frequency power spectrum for spherical and plane waves 
are given by equations (32) through (35) . The values of f 
for 9 April, 20 April and 30 April, based upon respective 
average wind speeds of 1 m/sec, 1.5 m/sec and 4 m/sec are 
30 Hz, 41 Hz and 53 Hz, respectively. For comparison with 
the theoretical model (Figure 3) , typical spectra for these 
experiments were plotted vs. "Log ft" and are shown in 
Figure 37. The vertical axis is the percentage of power 
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at a particular frequency compared to the maximum power 
attainable as determined by the spectrum analyzer’s 
"Y-axis” voltage setting. 

There appear to be two significant differences 

between the empirical results of Figure 37 and Clifford's 

theoretical spectra of Figure 3. First, Clifford indicates 

that the rapid downward slope of the spectrum begins at 

log ft ~-0.1 for plane waves and log ft ~ 0.2 for spherical 

waves. Empirical results show the drop occurring at 

log ft ~ 3.0 for the (7-9) April and 20 April experiments. 

The spectrum of (29-30) April shows an overall negative 

slope for all frequencies. Secondly, the slope appears to 

-8/3 

decrease less rapidly than the theoretical (f) ' . The 

slope for (7-9) April is of the order of (f)"^/H and for 
20 April and (29-30) April, (f)"2/13. p ux ther data 
collection is necessary to examine if, and to what extent, 
various turbulence parameters affect the magnitude of the 
negative slope. It does appear, however, that Clifford's 
(f) t00 rapid a decrease. 

4. Differences in Magnitude of the Aperture Averaging 
Effect 

Examination of Figures 30 through 33 indicates 
that there are significant differences between aperture 
averaging effects on the relative power spectrum as seen 
under the various experimental conditions. In an attempt 
to determine a factor common to the variances and to the 
power reductions which were not directly proportional to 
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aperture diameter, parameters such as wind speed, , 

C N , and G(D) were examined. No factor was found 

that individually accounts for the variations. A detailed 
analysis of all factors concurrently was not undertaken 
due to time limitations. 
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V. SUMMARY AND CONCLUSIONS 



A. APERTURE AVERAGING EFFECTS ON SCINTILLATION 

The variance of log-amplitude fluctuations that would 

theoretically be observed if using a vanishingly small 

receiver can best be approximated by an aperture of diameter 
1/2 

(Z/k) . A diameter of this size eliminates aperture 

averaging effects present when using apertures larger than 
1/2 

(Z/k) , and that may be present when using aperture 

1/2 

diameters smaller than (Z/k) . Experiments conducted 

in an enclosed laboratory in inhomogeneous turbulence and 
field experiments conducted in homogeneous turbulence, 
both over land and in the marine boundary layer, support 
this prediction from the theoretical model of Lutomirski 
and Yura [Ref. 1] . 

Fried's theoretical model of aperture averaging 

[Ref. 21] predicts that 0, the aperture averaging factor, 

increases continuously as aperture size is decreased from 

an arbitrarily large aperture diameter. Empirical results 

indicate that the aperture averaging factor does not increase 

continuously as aperture size is reduced to a vanishingly 

small diameter, but increases continuously until an aperture 
1/2 

diameter (Z/k) is reached. A further reduction in 

1/2 

aperture diameter [below (Z/k) ] causes 0 to (1) reverse 

direction and decrease or (2) remain near the same value 
of 0 . 
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In general, the value of the aperture averaging factor 

increases as <0 £j > increases. Put another way, for an 

1/2 

aperture of diameter larger than (Z/k) ' , the aperture 

averaging correction factor decreases as the amount of 

turbulence increases. Under weak turbulence conditions, 

the effect of aperture averaging is quite significant. For 

strong scintillation conditions, the aperture averaging 

effect is decreased and becomes less significant. 

The ratio G (D ’ ) / [G (D) ] m j_ n gives the factor by which the 

measured signal variance (or other turbulence parameter) , 

1/2 

measured with an aperture of diameter D 1 , D* > (Z/k) , 

must be multiplied to represent the value that would have 

been measured had a vanishingly small receiver been used. 

It is important to note that the value of [^(D)] is a 

function of the largest aperture diameter used (D ) , and 

should not be viewed as an absolute value for a particular 

set of turbulence parameters. In general, the value of 

[G(D)1 . decreases as D„ is increased. 

L v ' J mm o 

B. APERTURE AVERAGING EFFECTS ON THE TEMPORAL -FREQUENCY 
POWER SPECTRUM 

An increase in aperture diameter appears to decrease 
the magnitude of power in the frequency range 
.5 KHz < f < 2 KHz of the temporal -frequency power 
spectrum for scintillation measurements. The difference 
in power at a particular frequency for two given aperture 
diameters differed significantly for different experimental 
and turbulence conditions. The reason for the disparity 
was not determined due to time considerations. 
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Spectra obtained from experiments in the corridor 

laboratory, over land, and in the marine boundary surface 

layer indicate a negative sloped decrease in power at 

higher frequencies for the temporal- frequency power 

-8/3 

spectrum less rapid than the (f) ' predicted by 

Clifford [Ref. 19] . Further investigation is required 
to determine the effect changes in certain turbulence 
parameters have on the form of the power spectrum. 

C. EXTENSION OF PROJECT 

Comparisons of [G(D)] m i n with <a iil > indicated that, 
for a given aperture size larger than that which approxi- 
mates a vanishingly small receiver, the aperture averaging 
correction factor is not a constant, but is dependent upon 
the level of turbulence present. Further experimentation 
is needed to compare [G(D)] m ^ values, for a given maximum 

reference diameter D , with <o nT > values, to determine the 

o Jti 

relationship between the level of turbulence and the aperture 
averaging correction factor for a given aperture diameter D q . 

It seems clear from this report and other published 
data that aperture averaging reduces the high frequency 
power content in the temporal -frequency power spectrum 
and eliminates certain high frequencies proportional to 
increasing aperture diameter. However, experiments conducted 
under different turbulence conditions and over different 
propagation paths indicate that the power spectrum assumes 
different general shapes. This needs further investigation 
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through the examination of power spectra obtained over the 
same propagation path under various turbulence conditions, 
and under the same turbulence conditions over various 
propagation paths. In conjunction with this, Clifford's 
theoretical model of the temporal -frequency power spectrum 
for plane and spherical waves should be examined for validity 
more thoroughly through empirical study. 



77 



SUMMARY MATRIX OF EXPERIMENTS CONDUCTED 
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Path 1 (Apt Ave) Aperture 

Date Location Length (m) Detectors 2(W/Freq Anal) Diameters (mm) 
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Illustration of Kolmogorov Turbulence Model 
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Kolmogorov Turbulence Power Spectrum 
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Figure 3 




log ft 



Clifford’s theoretical temporal power spectrum of log- 
amplitude fluctuations for plane- and spherical-wave 
propagation. 
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Figure 4 
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PIN-5DP Photodiode Detector Circuit 
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Figure 6 

Field Lens System for Avalanche Photodiode Detector 
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Figure 7 

Field Lens System for PIN-5DP Photodiode Detector 
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Figure 8 

Electronic Components: Aperture Averaging Experiment 
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Figure 9 

Corridor Experiment Log- Intensity Distribution 



a 

a 




88 



R / v Figure 10 

ACANIA Propagation Paths over Monterey Bay 
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Spanagel Hall (NPS) 



G(D) vs. Aperture Diameter D 
Corridor Run, (7-8) February 1976 
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Aperture Diameter D(mm) 



Figure 12 

G(D) vs. Aperture Diameter D 
Corridor Run, 10 February 1976 



O 




O 

O 



O 

00 



o o o o 

'O (NJ o 



CD 



91 



345 67 12 17 22.6 (30) (40) 

Aperture Diameter D(mm) D = 22.6 mm 



Figure 13 

G(D) vs. Aperture Diameter D 
Corridor Run (7-9) April 1976 
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Aperture Diameter D(mm) 



Figure 14 

Comparison of Corridor Runs of 
February, 10 February and (7-9) April 
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Aperture Diameter D (mm) 



G(D) vs. D/(A.Z)^ 2 : Tatarski Model 
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Figure 16 

G (D) vs. D/(4Z/k) 1/2 : Fried Model 
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D/(4Z/k)^ 2 : Fried Model 



Figure 17 

1/7 

G(D) vs. D/(Z/k) ' : Lutomirski and Yura Model 
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D/(Z/k) 



Figure 18 

G(D) vs. D/(Z/k) ly/2 
Navy Annex Run, 20 April 1976 
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“/(Z/W D MO- 



Figure 19 

G CD) vs. D/(Z/k) 1/2 

Corridor, 7 Feb and (7-9) Apr; Navy Annex, 20 Apr 
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Aperture Diameter D(mm) D Q = 24.3 mm 

D . = 2 mm 



G(D) vs. Aperture Diameter, D 

Monterey Holiday Inn to Hopkins Marine Station, 11 March 1976 
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Aperture Diameter D(mm) D Q = 40 mm 

D . = 14 mm 



Figure 22 

G(D) vs. Aperture Diameter 
R/V ACANIA, 30 March 1976 
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Aperture Diameter D (mm) D q = 38 mm 

D . = 10 mm 



G(D) vs. Aperture Diameter 
R/V ACANIA, 2 April 1976 
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Aperture Diameter D (mm) D Q - 40 mm 

D . - 10 mm 



Figure 24 

G (D) vs . D/ (Z/k) 1 / 2 

R/V ACANIA, 2 April, and Holiday to Hopkins, 11 March 1976 
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D/(Z/k) 



Figure 25 

G(D) vs. D/(Z/k) 1/2 
R/V ACANIA, 27 April 1976 
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D/ (Z/k) D = 40 mm 



Figure 26 

G(D) vs. D/(Z/k) 1/2 
R/V ACANIA, 29 April 1976 
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13 mm 



Figure 27 

G (D) vs. D/(Z/k) 1//2 
R/V ACANIA, 30 April 1976 
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D/ (Z/k) D = 40 mm 



Figure 28 

Empirical and Theoretical "0 vs. r-7-7 

(4Z/k) i/z 

C4Z/k) i/2 




^ 7 Feb (Corridor) q (7-9) Apr (Corridor) 
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Figure 29 



”0 vs. rrr" : ACANIA Experiments 
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Relative Power Spectra 
(7-9) April 1976, Corridor Experiment 
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Frequency (Hz) 



Figure 31 

Relative Power Spectra 
20 April 1976, Over Land Experiment 
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Figure 32 

Relative Power Spectra 
20 April 1976, Over Land Experiment 
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Figure 33 

Relative Power Spectra 
(29-30) April 1976, R/V ACANIA Experiment 
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Figure 34 

Typical Power Spectrum 
9 April 1976, Corridor Experiment 
Aperture Diameter = 24.3 mm 
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Figure 35 

Typical Power Spectrum 
20 April 1976, Over Land Experiment 
Aperture Diameter = 24.3 mm 

Peak frequencies recorded from spectrum analyzer digital readout. 
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Figure 36 

Typical Power Spectrum 
(29-30) April 1976, R/V ACANIA Experiment 
Aperture Diameter = 25 mm 
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Figure 37 

vs. % of Maximum dK Scale Setting 
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q 20 Apr zk (29-30) Apr (Dashed line: Recorder linearity limit) 



APPENDIX B 



TEST PROCEDURE AND PERFORMANCE CHARACTERISTICS 
OF THE TEXAS INSTRUMENTS TYPE TIXL69 LARGE-AREA 
SILICON AVALANCHE PHOTODIODE 

A series of tests was conducted to evaluate the 
detection capabilities and performance characteristics of 
the TIXL69 avalanche photodiode detector and to compare 
them with that of the PIN-5DP photodiode detector. The 
two detectors were evaluated at a signal intensity compara- 
ble to that anticipated during a typical at-sea data run, 

- 9 2 

i.e., about 1.4X10 W/cm (the intensity of a 5 mil He-Ne 

O 

(6328A) laser source at a distance of about .8 nautical 
miles and with a beam divergence of 7 milliradians) . 

The circuit used with the avalanche detector is 
illustrated in Figure Bl(a). The reverse bias voltage 
source was a Hewlett Packard Harrison 6515A DC Power 
Supply. 

Dark current as a function of reverse bias voltage 
was measured first to see if it compared with the manu- 

-9 -7 

facturer’s stated values of 3.5X10 A minimum to 1.1X10 A 
maximum. The test circuit used is shown in Figure Bl(b). 
During this experiment the avalanche detector was covered 
by black cloth to keep light from hitting the detector 
surface. There were no windows in the enclosed room, all 
room lights were out, and the room was completely dark. 

A flashlight was used for recording data. To ensure that 
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the voltage applied was actually occurring as a voltage 
drop across the avalanche photodiode, a Hewlett Packard 
419A Null Voltmeter was connected to the circuit as shown. 
The voltmeter indicated that the avalanche photodiode was 
indeed receiving the full potential drop. Using a Keithley 
Instruments 600B Electrometer to measure the current, it 
appeared that values in the nanoamp region were occurring. 
However, the current indicator was not stable and repetitive 
current readings at the same reverse bias voltage fluctuated 
greatly. It was decided that a more accurate ammeter was 
needed and a Keithley Instruments 410 Micro -Microammeter was 
obtained. Steady, repeatable values of current were now 
attained for specific reverse bias voltages. 

As shown in Figure B2, dark current was linear with 
applied reverse bias voltage from 20 volts through 157.5 

3 

volts, then increased by a factor of about 10 at an 
applied voltage of 160 volts. Dark current values were 
higher than that stated by the manufacturer by a factor of 
10 at low applied reverse bias voltages, decreasing to a 
factor of 10 higher at high values. For this particular 
detector, breakdown occurs at a reverse bias voltage of 
160 volts vice the manufacturer’s indicated 166 volts. 

The HP Harrison 6515A DC voltage source was then re- 
placed by a battery box consisting of a series of seven 
22.5 volt and one 1.5 volt batteries to be used with the 
avalanche detector during actual aperture averaging 
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experiments. The "Dark Current vs. Reverse Bias Voltage" 
experiment was repeated. The same equipment and circuit 
as described above was used and dark current values obtained 
with the battery box were identical to those obtained with 
the HP Harrison 6515A DC source. 

With the breakdown voltage known, the next test was to 
examine the "Signal- to-Noise Ratio vs. Reverse Bias Voltage" 
at signal intensities similar to those previously described 
under at-sea conditions. The experimental setup was as 
shown in Figure B3. The signal source was a . 5 mW He-Ne 

O 

(6328A) laser. The laser beam was passed through a +8 cm 
focal length lens followed by a signal chopper located just 
beyond the focal point of the lens, then proceeded a distance 
of 3.35 meters to the avalanche detector. The avalanche 
detector was mounted on a stand that permitted horizontal and 
vertical movement in order to maximize signal intensity 
appearing on the oscilloscope. The signal from the detector 
was sent to a PAR (Princeton Applied Research) Model 113 
Pre-Amp, then to a Monsanto AM- 6419/USM- 368 oscilloscope 
for display. External triggering was used on the oscillo- 
scope. To obtain the desired signal intensity, a Jodon 
Model PM-550 Digital Optical Power Monitor, placed at the 
location of the avalanche detector, measured signal intensity 
directly. With the uninterrupted laser beam incident on the 
Jodon detector at close range (about 10 cm) , the Power 
Monitor indicated a maximum power output of .39 mW. The 
diameter of the laser beam was 6 mm wide while the diameter 
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of the active detection area of the Jodon detector is 3 mm. 

Since the laser beam is theoretically (and appeared to be 

from visual inspection) Gaussian, roughly 87% of the output 

intensity was incident on the Jodon detector active area, 

covering it completely, and giving a source power output 

of 5.5X10 W/cm . At the position of the avalanche detector, 

the intensity had decreased, due to beam expansion from the 

- 3 2 

unfocused optics of the source output, to 1.0X10 W/cm . 

With the chopper added and running at normal operating 

- 4 2 

speed, the intensity reduced further to 5.7X10 W/cm . 

With the diverging lens in place, to complete the setup, 

intensity decreased to 3X10 W/cm . To reduce the intensity 

further, neutral density filters (Kodak Wrattan Filter (in 

B glass) N.D.4, N.D.3 and N.D.2) were used singly and in 

combination. With the N.D.4 filter inserted (reduces 

intensity by 10 ’^) between the laser and the +8 cm lens, 

- 9 2 

signal intensity was reduced to 3X10 W/cm , as desired. 

When using the PAR, it was found that a saturated 
signal occurred on the oscilloscope if PAR Gain setting 
was too high. A saturated signal was identifiable as a 
rectangular-shaped 8.75 volt signal with normally visible 
scintillation absent. 

On 26 January 1976 an experiment comparing "S/N vs. 
Reverse Bias Voltage" at a given signal intensity for the 
avalanche detector was conducted followed by a comparison 
of "S/N of Avalanche Detector vs. S/N of PIN-5DP Detector" 
at varying intensities. In both experiments the equipment 
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setup of Figure B3 was used. PAR settings were: L-F 

ROLL OFF- -100 , H-F ROLL OFF--300K, Gain- -Variable . The 
Harrison 6515A DC Power Supply was used to supply reverse 
bias voltage to the avalanche detector and a 6328A filter 
was used in front of the avalanche and PIN photodiode 
detectors. No lens system was used to focus the incoming 
light onto the detector active area for either detector. 
Table BI gives data collected during the M S/N vs. Reverse 
Bias Voltage" for the avalanche detector, shown graphically 
in Figure B4. In Table BI, note (1) the lower PAR Gain 
settings required to avoid saturation as V^ increases 
beyond 120 volts, and (2) the increase in noise voltage 
as V, . increases. From Figure B4 it can be seen that S/N 

D 13.S 

increases fairly linearly from a reverse bias voltage of 
80 volts through 159 volts. Table BII is a comparison of 
the performance of the PIN photodiode detector and the 
avalanche photodiode detector. Both detectors were located 
at the same distance from the signal source and both used 
the same PAR and oscilloscope for signal output. It is 
evident that the avalanche detector, when set at its 
optimum reverse bias voltage, gives a significant improve- 
ment in performance over the PIN photodiode detector, 
especially at the weaker intensities characteristic of 
actual field experiments. 

A responsivity curve was plotted to compare experimental 
results with the manufacturer's stated performance charac- 
teristics. From Figure B5 it can be seen that experimental 
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results follow very closely the predicted values with the 
following exceptions: (1) actual signal gains corresponding 

to manufacturer’s highest predictions were not attained and 
(2) the maximum reverse bias voltage that could be applied 
was less than manufacturer indications. 

In summary, the TIXL69 avalanche detector proved to 
give a significant improvement over the PIN-5DP photo- 
diode detector for low intensity signal detection, as 
anticipated. The maximum (and optimum), reverse bias voltage 
that can be applied is 159 volts, slightly less than previous- 
ly expected. Measurements of responsivity agree with 
manufacturer's stated values but do not attain the highest 
values predicted. 
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Figure B1 (a) 

Circuit Diagram: Avalanche Photodiode Detector 
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Figure Bl(b) 

Test Setup: Dark Current Experiment 
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Figure B2 

Dark Current vs. Reverse Bias Voltage 
21 January 1976 
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Reverse Bias Voltage (Volts) 



Figure B3 

S/N Experimental Setup for Comparing Avalanche 
and PIN-5DP Photodiode Performance 
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TABLE BI 



S /N VS. V BIAS 

'FOR AVALANCHE DETECTOR 

O 

26 January 1976, Source: .5mW, 6328 A He-Ne laser. 

Using chopped signal and diverging lens, as shown in 

- Q 2 

Figure B3, signal intensity using N.D.F 4 ~ 3X10 W/cm . 



V BIAS 


PAR 

GAIN 


V (NOISE) 


V (SIGNAL) 


S/N 


20V 


10K 


. 55V 


1.4V 


2.5 


50V 


10K 


.60V 


1.55V 


2.6 


80V 


10K 


.65V 


1.9 V 


2.9 


100V 


10K 


. 70V 


2.75V 
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120V 


10K 


1.0 5V 
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150V 


2K 
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155V 


IK 


1. oov 


6.0 V 


6.0 


159V 


200 


.70V 


4.35V 
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160V 




NO 
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26 January 1976 
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Reverse Bias Voltage (Volts) 



COMPARISON OF AVALANCHE AND PIN-5DP PHOTODIODE DETECTORS 
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Responsivity (yA/yW) 



Figure B5 



Comparison of Manufacturer's Predicted and 
Empirical Responsivities 

Solid Line = Manufacturer's ©= Empirical 
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APPENDIX C 



EQUIPMENT PERFORMANCE CHARACTERISTICS 

1. Precision Instruments, Model No. PI-6204, four-channel 
recorder : 

Input sensitivity, .05 to 10 volts rms ; 

Frequency response, DC to 10 KHz +3 dB at 37.5 ips 
on FM Record; 

Signal to Noise Ratio, 35 dB at 37.5 ips, 300 Hz 
to 100 KHz. 

2. Sanborn HP-3914, 14-channel Analog Tape Recorder: 

Maximum Input, 5 volts (peak to peak) ; 

Gain, unity from 0 to 850 Hz, -.23 dB at 1 KHz; 
Tape Speed 7.5 ips. 
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